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Abstract 



The spin transfer effect in ferromagnet-quantum dot (insulator)-ferromagnet Aharonov-Bohm (AB) ring 
system with Rashba spin-orbit (SO) interactions is investigated by means of Keldysh nonequilibrium Green 



function method. It is found that both the magnitude and direction of the spin transfer torque (STT) acting 



on the right ferromagnet electrode can be effectively controlled by changing the magnetic flux threading the 
AB ring or the gate voltage on the quantum dot. The STT can be greatly augmented by matching a proper 



magnetic flux and an SO interaction at a cost of low electrical current. The STT, electrical current, and spin 
current are uncovered to oscillate with the magnetic flux. The present results are expected to be useful for 



information storage in nanospintronics. 
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I. INTRODUCTION 



The spin transfer effect (STE) states that when the spin-polarized electrons flow from one 
ferromagnet (FM) layer into another FM layer with magnetization aligned by a relative angle, they 
may transfer transverse spin angular momenta to the local spins of the second FM layer, thereby 
exerting a torque on the magnetic moments that is usually coined as the spin transfer torque (STT). 
This important phenomenon was predicted independently by Berger and Slonczeski Jl],Q] in 1996 
and soon confirmed by experiments. Because the STE can be utilized to switch the magnetic 
state of the free FM layer in a magnetic tunneling junction (MTJ) or a spin valve by applying 
an electrical current instead of a magnetic field, it may be even more useful in writing heads 
for magnetic random access memory (MRAM) or hard disk drivers than the conventional tunnel 
magnetoresistance (TMR) and giant magnetoresistance (GMR) effects. In view of the potentially 
wide applications in nanospintronic devices, a number of works on the STE have been done for 
different systems both theoretically and experimentally f 3I— 161. 

On the other hand, the quantum dot (QD) has received much attention in the past decades, and 



a lot of advances have been made in this particular field (e.g. Refs. ill 



2 IIP . For a semiconductor 



QD, as the spin-orbit (SO) interaction is usually not negligible, some interesting phenomena re- 
lated to the SO interactions, such as the bias-controllable intrinsic spin polarization in a QD [12211 
and the interplay of Fano and Rashba effect B23I1 . can be observed. Almost twenty years ago, Datta 
and Das predicted a spin transistor based on the Rashba SO 112411 . showing that the SO interactions 
may be important in the semiconductor spintronics. However, the effect of the Rashba SO interac- 
tion on the STE is still sparsely studied. In this paper, we shall take the FM-QD (insulator, I)-FM 
Aharonov-Bohm (AB) ring system as an example to investigate how both the SO interaction and 
the magnetic flux affect the spin-dependent properties of the system by means of the nonequilib- 
rium Green function method. We have found that the magnitude and direction of the STT can be 
easily controlled by changing the gate voltage V g on the QD or the magnetic flux <p through the 
ring if both the SO and electron-electron (e-e) interactions in the QD are considered, which might 
be useful in information storage. 

The other parts of this paper are organized as follows. In Sec. II, a model is proposed, and 
the relevant Green functions are obtained in terms of the nonequilibrium Green function method. 
In Sec. Ill, the spin-dependent properties of STT in the system under interest are numerically 
investigated, and some discussions are presented. Finally, a brief summary is given in Sec. IV. 
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II. MODEL AND METHOD 



The system under interest is depicted in Fig. 1. Two FM leads spreading along the z axis are 
weakly coupled to an insulating (I) barrier and a semiconducting QD, forming an AB ring. The 
left (L) FM electrode with the magnetization along the z axis is applied by a bias voltage —V/2, 
while the right (R) electrode with the magnetization along the z' axis that deviates by an angle 6 
from the z axis is applied by a bias voltage V/2. Assume that the QD is made of a two-dimensional 
electron gas in which the electrons are strongly confined in the y direction by a potential V(y). 
Due to dV/dy 3> dV/dx and dV/dz, we have W(y) ~ y(dV/dy), where y is the unit vector 
along the y axis. If V(y) is asymmetric to y = 0, both Rashba SO and e-e interactions on the QD 
should be considered. Since the electronic transport of the device along the z axis is much more 
dominant than that along other two dimensions, the device under interest can be treated as a quasi 
one-dimensional system. Sun et al. y2_5J] have carefully analyzed the SO Rashba interaction and 
found that (i) the Rashba SO interaction can be separated into two parts, H Rl and H R2 , namely 

y 

H so = —■ [a(x)(a xp) + (?x p)a(x)\ = H Rl + Hr 2 , (1) 
H Rl = —[a(x)a z p x + a z p x a(x)}, (2) 

h R2 = — L -J-^ i ; (3) 

(ii) by choosing a suitable unitary transformation, H Rl can give rise to a spin-dependent phase 
factor in the tunneling matrix element between the leads and the QD, while Eq. (3) can be written 



in the second-quantization form as [25]: H R2 = XX^mn^m^nt + h.c), that causes a spin-flip 

mn 

term with strength t™ n in the QD, where m and n are quantum numbers for the eigenstates of 
electrons in QD; (iii) since the time-reversal invariance is maintained by the Rashba SO interaction, 
tmn = ~t n °m m ^ ^tin = 0' which suggests that the spin-flip scatterings only occur between different 
levels in the QD. In the present work, for simplicity, we shall consider the case with a single-level 



QD as in some previous works [|23Ll25|l . where no interlevel spin-flip scattering happens in the QD. 
Thus, H R2 equals to zero. Suppose that a(x) is independent of the coordinates in the scattering 
region, and a magnetic flux penetrates into the AB ring. The Hamiltonian of the present system is 
given by 

H = Hqd + H/3 + H T , (4) 
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Hqd = ^2 edd t d ° + u n t n h ( 5 ) 



Hp — y^g/3fcgQgfca Q /3fc^^ (6) 

0k,a 



k,cr 



+ Y^[ t Ldat ka d a + h.c.} 

k,a 

+ Yl^ tLR ( C0S -j; a+ Rkv ~ ° Sil1 -^ a Rka) a Lka + h.C.}, (7) 



k.a 



where ap ka and d a are annihilation operators of electrons with momentum k and spin a (=f, |) 
in the /3 (= L, R) electrode and in the QD, respectively, ep k(y = e k + aMp — eVp is the single- 
electron energy for the wave vector k with the molecular field Mp in the electrode (3, Ed is the 
single-electron energy in the QD, U represents the on-site Coulomb interaction between electrons 
in the QD, tpa is the tunneling matrix element of electrons between the j3 electrode and the QD, 
t LR is the tunneling matrix element of electrons between L and R electrodes through the insulating 
barrier, n a = c+c a , and 7 = k R d with k R = am*/h 2 , a = (&(y) \{d/dy)V{y)\ ^{y)), m* the 
effective mass of electrons and d the thickness of the middle region. The magnetic flux $ threading 
the AB ring is related to the phase factor by 4> = 27r$ / <3> , where $ is the flux quantum. It should 
be noted that the magnetic flux threading the AB ring generally includes two contributions, one 
generated by the FM leads that may be small and constant, and the other from the external magnetic 
field that can be varied to adjust the phase factor <p. 

The transverse component of the total spin in the right FM lead can be written as yQ 

& = 2 (^fct' a Rki) 

k 

= 2 ^( a Rkt a Rki + a Rki a m)> (8) 
k 

where S is written in the x'y'z' coordinate frame. The spin torque, namely, the time evolution rate 




of the transverse component of the tota 



( [H, S]). According to Refs. J3,|26|,|27j], the right FM layer gains two types of torques: one is the 



spin of the right FM lead, can be obtained by dS/ dt 
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FIG. 1 : A schematic layout of FM-QD(I)-FM AB ring system. The electrons flow towards the x axis. 



equilibrium torque caused by the spin-dependent potential, and another is from the tunneling of 
electrons that is in what we are interested. After cautiously separating the current-induced torque 
from the equilibrium one, the STT is given by 











r = -Ke{t M e*V^cos-G*< (t, t) - t M e"^e^ sin -G< (t,t) 



+t RL cos -G^ m (t, t) - t RL sin ^Gf m (t, t) 




-( 

2 











tjufi^e* cos -G<(t, t) + tuffe* sin -G<(t, t) 











+t RL cos -Gi W (t, t) + t RL sin -GU(i, t)}. 



(9) 



From Eq. (9), it is clear that the current-induced STT can be obtained as long as we get the 
lesser Green functions G < . In what follows we shall use Keldysh's nonequilibrium Green function 
technique to determine all lesser Green functions J28]. These functions are closely related to the 
retarded Green functions defined by 



k' k 

G^(t,t>) = -i9(t-t')({J2^ka(t),d+(t')}), 

k 

= -<e(t-0(K(*),^(0}>, 

where {A, B} denotes the anticommutation relations, and (A) stands for the thermal average. By 
using the equation of motion, the retarded Green functions can be obtained by Dyson equation 
G r = g r + g r T l r G r , where g r is the retarded Green function for decoupled systems, and S r is the 
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self-energy of electrons. To obtain G r , the decoupling approximations similar to those in Refs. 



[5, 



28 



3 lh for the equations of motion of Green functions should be made. As the associated 



equations for Green functions are quite lengthy, we shall not repeat them here for conciseness. 
The lesser Green function G < can be calculated straightforwardly from the Keldysh equation 



G< = (l + G r Z r )g < {l + Z a G a ) + G r Y, < G a 
= G r g r - 1 g < g a - 1 G a + G r lu < G a . 



In the present case, = 0, and g r 1 g < g a 1 is diagonal 



(10) 



9 r - 1 9 < 9 a - 1 



2ifL(e)/npLi 



2if R (e)/irp Rr 



\ 



2if R (e)/7ip Rl 



\ 



The electrical current is given by 



(11) 



I = h + h, 



J t = -^ { t M e-^cos^< w (M )+ t M e^s m ^< 
+t R L cos -G< 



Q 

e bV cos -G^ Rt {t,t) +t Rd e 

9 9 
g^it-Kt^M) +tRL sin -G^ m {t,t)}, 



h = -^e{t Rd e^e^ C os 9 -Gl Ri (t,t)-t Rd e- l ''e^sm d - 
9 



(12) 



(13) 



The 



sm-G^it^) 

9 9 
+t RL cos -G^ m (t,t) - t RL sin -G^ m (t,t)}. (14) 

spin current is defined by a difference between the electrical currents of spin up and down, 

T — T t nz\ 



h=h-h- 



To get the physical quantities of interest, the above-mentioned equations will be solved numeri- 
cally in a self-consistent manner. 
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III. RESULTS AND DISCUSSIONS 



It has been shown that when the incident electrical current is larger than a critical value, the STT 
can switch the direction of the magnetization of the free FM layer clockwise or anticlockwise de- 
pending on the direction of the incident electrical current II29I-I3 ill . In the present case, the positive 
STT tends to push the spins in the right FM electrode aligning antiparallel with the magnetization 
of the left FM electrode, while the negative STT may cause a reverse orientation of the magnetiza- 
tion in the free FM layer. In order to properly incorporate the STE into a functionalized spintronic 
device, both the direction and magnitude of the STT should be taken into account. For simplicity, 
in the following parts we will assume that in most cases the left and right FM electrodes have the 
same spin polarization P L = P R = P = 0.5, and the angle 9 between z and z' axes is tt/3 through- 
out the paper unless specified. We take 7 = ^ and T = T L ^(P = 0) = T L ^ R )^(P = 0) 
as scales for the electrical and spin currents as well as the STT and energy, respectively, where 



Tccr^) = 2tt ^2 ka \t a d\ 2 S(e — £fc„). In accordance with Refs. 
SO interaction constant is a ~ 3 x lO~ n eVm, k R = m*a/h 2 ps 0.015/nm for m* = 0.036m e , 
the typical length of QD is 100 ran, U = 5T , and 7 can be n/2 or larger. 

As the STE only exists in the noncollinear case, in contrast to previous works where the electri- 
cal and spin currents were discussed only in collinear cases (9 = or it) (e.g. |23i|250), let us first 



5. E- 



3411 . we assume that the Rashba 



look at the angular dependences of the STT, the electrical current and spin current for different 
magnetic flux <p and Rashba SO interaction 7. The results are given in Fig. 2, where 7 = it/ 2 and 
Ed = 1 in Figs. 2(a)-(c), and = n/2 and e d = 1 in Figs. 2(d)-(f). It can be observed that the STT 
has a sine-like relationship with the relative angle 8, and the direction and magnitude of the STT 
are clearly influenced by both the magnetic flux <p and Rashba SO interaction 7. In the absence of 
either <\> or 7, the STT remains negative (anticlockwise). For the simultaneous presence of </> and 
7 (greater than it /A), the STT becomes positive (clockwise). This fact reminds us that we may 
apply the magnetic flux to change the direction of the STT, thereby being capable of manipulating 
the magnetic state of the free FM layer, which might be useful for information storage and for 
designing the memory element. The electrical current I decreases with increasing 9 in the absence 
of either <p or 7, indicating a spin-valve effect, while it increases with 9 in the presence of both <p 
and 7 (greater than it /A), giving an anti-spin- valve effect [e.g. or 7 = tt/2 in Figs. 2(b) and (e)]. 
This property differs obviously from the conventional FM-I-FM or FM-QD-FM systems without 
considering the SO interactions where the electrical current always decreases with increasing 9. 
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FIG. 2: (Color Online) The dependence of (a) the spin transfer torque r, (b) the electrical current /, and 
(c) the spin current I s for different magnetic flux <ft. The 9 dependence of (d) the spin transfer torque r, (e) 
the electrical current /, and (f) the spin current I s for different Rashba SO interaction 7. 

The spin current shows a feature similar to the electrical current. From these calculated results 
presented in Fig. 2, we can find that for a given Rashba SO interaction 7 (magnetic flux 0), the 
angular dependent STT, electrical current and spin current exhibit distinct behaviors for different 
magnetic flux (Rashba SO interaction). 

Figure 3 shows the magnetic flux </> dependence of the STT, electrical current and spin current 
for different Rashba SO interactions. We can see that with increasing 0, the STT, electrical current 
and spin current oscillate differently for various Rashba SO interactions 7. The larger the SO 
interaction 7 is, the more complex the oscillations are. For the QD energy level e<i = — 1 and 
7 = 7r/2, the STT shows a maximum around = 3n/2, while the electrical current exhibits 
minima around the same 0. Therefore, we may be able to use a lower current to change the 
magnetic state of the free FM by adjusting the magnetic flux penetrating into the AB ring. It is 
favorable for the spintronic devices, because a larger current may cause more heating, while the 
heating should be reduced as small as possible for better functions of the device. In addition, it 
can be found that the STT is closely related to the spin current, as the dips and peaks of Figs. 3(a) 
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FIG. 3: (Color Online) The magnetic flux dependence of (a) spin transfer torque r, (b) electrical current /, 
and (c) spin current I s for different Rashba SO interaction 7, where Ed = —1, and = tt/3. 

and (c) appear almost at the same positions. 

The energy level Ed of electrons in the QD, that offers resonant tunneling channels for spin- 
polarized electrons from the left FM electrode to the right FM one, has also effects on the magnetic 
flux 4> dependence of the STT, electrical current and spin current. The results are presented in Fig. 
4 for 7 = 7r/2. It is unclosed that for different r, / and I s exhibit different features, and 
oscillate with <\> in general. When Ed = 0, r, /, and I s are mirror symmetrical to (ft = n, and r 
is always negative. For positive and negative E d , r and I s have just opposite properties: the peaks 
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FIG. 4: (Color Online) The magnetic flux dependence of (a) the spin transfer torque r, (b) the electrical 
current /, and (c) the spin current I s for different energy levels e d , where = 7r/3, and 7 = ir/2. 

at = tc/2 (dips at = 3ir/2) for e d = 1 correspond to the dips (peaks) for e d = —1 at the 
same <fr, but the curves for positive and negative Ed intersect at = n, as shown in Figs. 4(a) and 
(c). It hints us that by changing the gate voltage that is usually utilized to alter the energy levels 
Ed in the QD, one can adjust the STT. For example, when <p = vr/2, if we change Ed from — 1 
to 1, the STT changes from 0.06 anticlockwise to 0.075 clockwise. As the gate voltage is easier 
than the magnetic flux to control, the present observation may offer a useful way to manipulate the 
magnetic state of the free FM layer. The electrical current also displays quite different oscillating 
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FIG. 5: (Color Online) The magnetic flux dependence of the spin transfer torque for different spin polariza- 
tions at 9 = 7r/3, and 7 = tt/2. 

behaviors for e^ — \ and —1, which is shown in Fig. 4(b). 

Why can the STT be controlled by changing the magnetic flux and gate voltage? Because the 
transmission probability of the spin-up electrons is proportional to cos(9 + + 7) and that of 



spin-down electrons is proportional to cos (9 + <p — 7) [|25|1 . The spin-up and spin-down electrons 
have different transmission probabilities if 7 is nonzero, leading to oscillations of the STT with 
magnetic flux </>. On the other hand, the STT is intimately related to the electrical current / [51], 
and the magnitude of I depends on the energy level of QD, so it is reasonable that the STT can 
be manipulated by adjusting the gate voltage. J-j. and I± give rise to opposite STT on the right FM 
layer. Since If and Zj. oscillate for various combination of <p and 7 in different ways, the STT may 
reach the maximum while I is in its minimum when <p and 7 take proper values. 

In addition, the spin polarization Pl,r of the FM electrodes has also effects on the magnetic 
flux dependence of the STT, as shown in Fig. 5. Generally, with increasing Pl,r, the STT shows 
qualitatively similar behaviors for positive and negative e d . When P L = P R = P, as shown in 
Figs. 5(a) and (b), the larger the polarization P, the smaller the peaks of the STT. For different P, 
t has obvious changes when < <p < 3tt/2 for Sd = 1, and when tt/2 < < n for Ed = — 1. 
When Pr and Pl are different, e.g. Pr = 0.5 and Pl = 0, 0.3, 0.7, the larger Pi is, the more 
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downward the curves move, as indicated in Figs. 5(c) and (d). It is interesting that even the left 
electrode becomes spin unpolarized (Pl = 0), the STT as a function of still behaves a sine- 
like curve and retains almost intact for different spin polarizations Pr [Figs. 5(e) and (f)]. The 
existence of the STT at P L = demonstrates that even if the left electrode is a normal metal 
(NM), the unpolarized electrons from the left NM lead flowing into the AB ring system with an 
QD encompassed by a magnetic flux can become spin-polarized before entering into the right 
FM electrode. It is apprehensible , because owing to the Rashba effect, the spin-up and spin-down 
electrons pass through the AB ring system at different transmission probabilities, as discussed 
above. When these spin-polarized electrons flow into the right FM layer, they may transfer some 
spin angular momenta to the local spins of the right FM electrode, thereby giving rise to the STT. 
In this case, if = 0, the STT becomes negligibly small. In the above analysis, we have presumed 
that the spin relaxation time of electrons is greater than that of the tunneling time. Thus, to ensure 
the feasibility of experimental observation, one must choose proper materials as FM electrodes 
and QD, and design a viable ring system to meet with the above requirements. It is interesting 
to note that a similar mesoscopic ring system was proposed, where some material parameters 
were discussed for possible experimental implementation (13511 that may be insightful for choosing 
proper materials for designing the present ring system. 

The effect of Rashba SO interaction 7 on the STT, electrical current and spin current is shown 
in Fig. 6 for different magnetic flux 0. With increasing 7, when = 0, the STT is always negative 
and goes down non-monotonously. When = 7r/4 or 7r/2, r goes up from negative to positive, 
reaches a round maximum, and then decreases, as depicted in Figs. 6(a). This result implies that 
the STT can be enhanced remarkably by matching with proper 7. The 7 dependences of the 
electrical current and spin current show different behaviors for various 0, as presented in Figs. 
6(b) and (c). With increasing 7, for = 0, both I and I s increase; for = 7r/4, / first decreases to 
a round minimum, and then goes up, while I s declines slowly; for = tt/2, the situation becomes 
reverse, i.e., / decreases dramatically, while I s first declines and then goes up. In a word, the 
Rashba SO interactions have various effects on the STT, electrical current and spin current. 

Finally, the bias voltage dependences of the STT, electrical current and spin current are studied 
for different 7 and 0, as shown in Fig. 7. In the simultaneous presence of 7 and 0, e.g. 7 = = 
7r/2, with increasing the voltage, the STT first increases almost linearly, reaches a peak, and then 
decreases slowly. After reaching zero, it starts to increase again in a different direction. In the 
absence of either 7 or or both, r is negative, and decreases non-monotonously with increasing 
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FIG. 6: (Color Online) The Rashba SO interaction 7 dependence of (a) the spin transfer torque r, (b) the 
electrical current I, and (c) the spin current I s for different magnetic flux <ft, where 9 = it/ 3, and = 1. 

the bias voltage, as indicated in Fig. 7(a). For various combinations of 7 and 0, the electrical 
current / exhibits qualitatively similar behaviors, which increases overall in a non-ohmic way 
with increasing the bias [Fig. 7(b)]. For 7 = c/>, I s remains almost constant at a small bias. When 
the bias passes a threshold, it increases linearly with the increase of V. For 7 (p, I s grows up 
almost linearly despite of small shoulders at a low bias, as displayed in Fig. 7(c). From Figs. 7(a) 
and (c), we can see that the shoulder structure of the STT and the threshold of the spin current 
appear around eV = 2ed, where the resonant tunneling happens. It is not surprising that the 
resonant tunneling has influences on the spin-dependent transport of the system. However, it is 
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FIG. 7: (Color Online) The bias voltage dependence of (a) the spin transfer torque r, (b) the electrical 
current /, and (c) the spin current I s for different 7 and <fi. 

more important when 7 = 0, while it is negligible when 7 7^ 0. 
IV. SUMMARY 

By means of the Keldysh nonequilibrium Green function method, we have investigated the 
STE in the FM-QD(I)-FM ring system with Rashba SO interactions. It has been found that both 
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the direction and magnitude of the STT are affected by the magnetic flux and the Rashba SO 
interactions. When the SO interaction is strong enough, the STT acting on the spins of the right 
FM electrode can be remarkably enhanced by matching the magnetic flux through the AB ring, 
which makes it is possible to readily manipulate the magnetic state of the free FM layer at a cost 
of lower electrical current. This property is quite expected for nanospintronic devices where the 
excessive heating generated by the electrical current should be avoided as much as possible. It has 
also been uncovered that by adjusting the gate voltage acting on the QD, both the magnitude and 
the direction the STT can be changed, which gives an alternative way to manipulate the magnetic 
state of the free FM layer. In addition, it is interesting to observe that the STT can also be increased 
by the magnetic flux through the ring or the gate voltage on the QD even if the left FM lead is 
changed to a NM. 

We would like to mention that the results presented in this paper provide useful information 
for designing practical spintronic devices based on the STE. Such a ring layout can be used either 
as a memory element with a low driving current or as a magnetometer to measure weak magnetic 
fields, because the tunnel current depends sensitively on the magnetic flux threaded the ring. On 
the other hand, the tunnel current or the magnetic state of the free FM layer are affected by the 
Rashba SO interaction, and one may inversely enable to estimate the magnitude of the Rashba SO 
interaction on the QD by means of such a ring apparatus. We expect that the present theoretical 
findings could be tested experimentally in future. 
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